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Abstract 

In two previous papers the author introduced a multiplication of distributions in one di- 
mension and he proved that two one-dimensional Dirac delta functions and their derivatives 
can be multiplied, at least under certain conditions. Here, mainly motivated by some en- 
gineering applications in the analysis of the structures, we propose a different definition of 
multiplication of distributions which can be easily extended to any spatial dimension. In 
particular we prove that with this new definition delta functions and their derivatives can 
still be multiplied. 



1 Introduction 



In the literature several examples of multiplication of distributions exist, more or less interesting 
and more or less useful for concrete applications. In [2] and [3j we have proposed our own 
definition of multiplication in one spatial dimension, d = 1, and we have proved that two 
or more delta functions can be multiplied and produce, as a result, a linear functional over 
T>(M). However, our definition does not admit a natural extension to d > 1. This is a strong 
limitation, both mathematically and for concrete applications: for instance, it is known to the 
civil engineers community that the problem of fracture mechanics may be analyzed considering 
beams showing discontinuities along the beam span, [HE]. The classical approach for finding 
solutions of beams with discontinuities consists in looking for continuous solutions between 
discontinuities and imposing continuity conditions at the fractures. In [HE] a different strategy 
has been successfully proposed, modeling the flexural stiffness and the slope discontinuities as 
delta functions. However, in this approach, the problem of multiplying two delta functions 
naturally arises, and this was solved using the definition given in [2]. This application proved 
to be useful not only to get a solution in a closed form but also to get numerical results in a 
reasonable simple way, when compared with the older existing approaches. These very promising 
results, however, have been discussed only in d = 1 since the mathematical framework, which 
is behind the concrete application, only existed in one dimension, [3]. It is not surprising, 
therefore, that extensions of our procedure to higher spatial dimensions has been strongly urged 
in order to consider the same problem for more general physical systems, i.e. for physical 
systems in higher dimensions like two or three-dimensional beams, in which the fractures can 
be schematized as two or three-dimensional delta functions. 

In a completely different field of research the same necessity appeared: for instance, in the 
analysis of stochastic processes the need for what is called a renormalization procedure for the 
powers of white noise, where the white noise is essentially a delta function, has produced many 
results, see pQ. Also, applications to electric circuits exist which are surely less mathematically 
oriented, [6], and again are based on the possibility of giving a meaning to 5 2 (x). With these 
motivations we have considered the problem of defining the multiplication of two distributions in 
more than one spatial dimension. However, our original definition cannot be easily extended to 
d > 1, because the regularization proposed in [7] and used in [21 [3] does not work without major 
changes in this case. For this reason we propose here a different definition of multiplication, 
which works perfectly in any spatial dimension. Moreover, this new approach returns results 
which are very close to those in [21 [3J, for d = 1. 

The paper is organized as follows: 

in the next Section we briefly recall the main definitions and results of [2] and [3j; 
in Section 3 we propose a different definition of multiplication in one dimension and we show 
that results which are not significantly different from those of Section 2 are recovered; 

in Section 4 we extend the definition to an arbitrary spatial dimension and prove that with 
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this new definition two delta functions can be multiplied. 

2 A brief resume of our previous results 

In [21 [3] we have introduced a (family of) multiplications of distributions making use of two 
different regularization discussed in the literature and adapted to our purposes. Here, for com- 
pleteness' sake, we briefly recall our strategy without going into too many details. 

The first ingredient was first introduced in |7j, where it was proven that, given a distribution 
T with compact support, the function 

T°(z)^^-T-(x-z)- 1 (2.1) 

exists and is holomorphic in z in the whole z-plane minus the support of T. Moreover, if T(x) is 
a continuous function with compact support, then T re( i(x, e) = T°(x + ie) — T°(x — it) converges 
uniformly to T(x) on the whole real axis for e — > + . Finally, if T is a distribution in D'(M) with 
compact support then T rec i(x, e) converges to T in the following sense 

/oo 
T red (x,e)*(x)dx 
-oo 

for every test function \P G D(M). 

As discussed in [7], this definition can be extended to a larger class of one-dimensional 
distributions with support not necessarily compact, V(M), while it is much harder to extend the 
same definition to more than one spatial dimension. 

The second ingredient is the so-called method of sequential completion, which is discussed 
for instance in [5], and it follows essentially from very well known results on the regularity of 
the convolution of distributions and test functions. Let <3? G P(R) be a given function with 
supp C [—1, 1] and J R <3?(x) dx = 1. We call 5— sequence the sequence S n , n G N, defined by 
S n ( x ) = n&(nx). Then, VT € D'(M), the convolution T n = T * 6 n is a C°°— function for any 
fixed n G N. The sequence {T n } converges to T in the topology of T>' , when n — > oo. Moreover, 
if T(x) is a continuous function with compact support, then T n (x) converges uniformly to T(x). 

We are now ready to recall our original definition: for any couple of distributions T, S G 
P'(R), Va,/3 > and G T>(M) we start defining the following quantity: 




where 

S^(x)^(S*S^)(x), (2.3) 

with 5^\x) = n /3< l ) (n /3 x), which is surely well defined for any choice of a,/?, T, S and ^. 
Moreover, if the limit of the sequence ((S ® T)^"'' 3 ^(^')| exists for all \l/(z) G T>(K), we define 
(S®T) (0hP) (V) as: 

(^®r)(a,/3)(*) = I™ (S®T)£»'fl(¥) (2.4) 
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Of course, as already remarked in [2], the definition (j2.4f) really defines many multiplications 
of distributions. In order to obtain one definite product we have to fix the positive values of 
a and /3 and the particular function $ which is used to construct the (^-sequence. Moreover, 
if T(x) and S(x) are two continuous functions with compact supports, and if a and /3 are any 
pair of positive real numbers, then: (i) Tn (x) S re d(x, \) converges uniformly to S(x) T{x); (ii) 

V*(x) G V(R) 4(T® S) (a ,fl(*) = IZo T ( x ) S ( x ) dx - 

It is furthermore very easy to see that the product (5 ®T)r a ^-\ is a linear functional on T>(M) 
due to the linearity of the integral and to the properties of the limit. The continuity of such a 
functional is, on the contrary, not obvious at all, but, as formulas (|2.6p - (|2.1ip show, is a free 
benefit of our procedure. 

We now recall few results obtained in [21 [3] . 
If we assume <3> to be of the form 

*(*) = { €'-P{^}> N<! (2 . 5) 
\ 0, \x\ > 1. 

where m is an even natural number and N m is a normalization constant which gives 3>(x) dx = 
1, and we put Aj = -^p- dx, whenever it exists, we have proved that: 
if m > 1 then 



±A 2 5, a = 2/3 
0, a > 2/3; 



if m > 2 then 



if m > 3 then 



Also, if m > 3 then 



If m > 4 then 



Finally, if m > 5 then 




(2.6) 



Va > 3/3; (2.7) 

=M 4 <5, a = 4/3 . , 

(2.8) 

(2.9) 
(2.10) 
(2.11) 
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a = 6/3 



0, a > 6/3. 

It is worth stressing that formula (|2.6|) . for a > 2/3, coincides with the result given by 
the neutrix product discussed by Zhi and Fisher, see [9]. Also, because of our technique which 
strongly relies on the Lebesgue dominated convergence theorem, LDCT, the above formulas only 
give sufficient conditions for the product between different distributions to exist. In other words: 
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formula (|2.6[) does not necessarily implies that (5 ® ^)( a ,/3) does not exist for a < 2/3. In this 
case, however, different techniques should be used to check the existence or the non-existence of 

More remarks on this approach can be found in [2] where, in particular the above results are 
extended to the product between two distributions like 5® and 5^ for generic I, k = 0, 1, 2, . . .. 
In [3] we have further discussed the extension of the definition of our multiplication to more 
distributions and to the case of non commuting distributions, which is relevant in quantum field 
theory. 



3 A different definition in d = 1 

The idea behind the definition we introduce in this section is very simple: since the regularization 
T — > T Te d cannot be easily generalized to higher spatial dimensions, d > 1, we use twice the 
convolution procedure in (12. 31) . T — > Tn =T * Sn , with different values of a as we will see. 

Let <I>(x) G T>(M) be a given non negative function, with support in [—1,1] and such that 
J K <I>(x) dx = 1. In the rest of this section, as in [21 [3] , we will essentially consider the following 
particular choice of ^(x), 

*(*) = ( €- 6Xp ^>' N<1 (3.1) 
{ 0, \x\>l. 

where m is some fixed even natural number and N m is a normalization constant fixed by the 
condition f_, $(x) dx = 1. As we have discussed in [2] the sequence 5n (x) = n a Q(n a x) is a 
delta- sequence for any choice of a > 0. This means that: (1) 5^(x) — > 5(x) in X>'(R) for any 
a > 0; (2) for any n S N and for all a > if T(x) is a continuous function with compact support 
then the convolution T^\x) = (T * 5^)(x) converges to T(x) uniformly for all a > 0; (3) if 
T(x) G D(R) then ri Q) (x) converges to T(x) in the topology of V(R); (4) if T e P'(M) then 
Tn a \x) is a C°° function and it converges to T in D'(M) as n diverges independently of a > 0. 

We remark here that all these results can be naturally extended to larger spatial dimensions, 
and this will be useful in the next section. 

Our next step is to replace definition (|2.4p with our alternative multiplication. To begin 
with, let us consider two distributions T,S G Z?'(R), and let us compute their convolutions 
Ti a) (x) = (T*5 [ n a) )(x) and S$P(x) = (T*5 { ®)(x) with ^ q) (x) = n a <S>(n a x), for a,f3 > to be 
fixed in the following. Tn°\x) and S^\x) are both C°° functions, so that for any ^(x) G D(M.) 
and for each fixed n G N, the following integrals surely exist: 

(SQTt^W = \( \s^(x)T^(x) + S^\x)T^{x)\ *(z)dx, (3.2) 
(5 d r)M(*)= / 4")(x)T^)(x)*(x)dx, (3.3) 



5 



(S© e *T)^)(*)= / sW(x)I< a Hx)*(x)dx. (3.4) 

JR 

The suffix d stands for direct, while ex stands for exchange. This is because the two related 
integrals remind us the direct and the exchange contributions in an energy Hartree-Fock compu- 
tation, typical of quantum many-body systems. It is clear that if S = T then the three integrals 
above coincide: (S S) { n' P) (V) = (S Q d S)^^) = (S Q ex S)^^), for all a,f3,^,n. In 
general, however, they are different and we will discuss an example in which they really produce 
different results when n — > oo. We say that the distributions 5 and T are 0-multiplicable, and 
we indicate with (S T)r a m their product, if the following limit exists for all ^(x) G D(R): 

(SQT) M (*) = lim(S0T)^)(*). (3.5) 

v ' n,oo 

Analogously we put, when the following limits exist, 

(S0 d T) (a>/9) (*) = lim (S0 d T)W)(f). (3.6) 

and 

(S Q ex T) {a>/3) (y) = lim (S Q ex r)^(*). (3.7) 

Again, it is clear that, whenever they exist, (S S) r a p\ ("if) = (S 5 , )( / 3 iQ )( x I') = (S Q d 
s )(a,f3){^) = (S Qex £) («,/?) W, for all ^(x) G V(R) while they are different, in general, if 
S ^ T. Of course, the existence of these limits in general will depend on the values of a and f3 
and on the particular choice of <J>(x). For this reason, as in [21 [3] , we are really defining a class 
of multiplications of distributions and not just a single one. We have already discussed in [2] a 
simple example which shows how these different multiplications may have a physical interpre- 
tation in a simple quantum mechanical system. We will discuss further physical applications of 
our procedure in a forthcoming paper. 

We now list a set of properties which follow directly from the definitions: 

1. if S(x) and T(x) are continuous functions with compact support then 

(S0T) W) (f) = (S0 d T) W) (*) = (5 e ,T) M (f) = / S(x)T(x)^(x)dx, (3.8) 

JR 

for all $(ar) G V(R) and for all a,/3 > 0. 

2. for all fixed n, for all a, (3 > and for all ^(x) G D(R) we have 

(5 Q d T)fr® (tf ) = (S Q ex T)^ (*) (3.9) 

3. for all n, a, /3 > and for all \&(x) G P(R) we have 

(5 Q d Tftf* (*) = (T d S)f q) (¥), and (5 e:c T)^ (tf) = (T Q ex S)^ a) (*) 

(3.10) 
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4. given S, T € T>'(M) such that the following quantities exist we have 

(s o r) (a>jB) (*) = ^{(^0d r) (a>jB) (*) + (5 e , r) (ai/J) (*)} (3.11) 

We will now discuss a few examples of these definitions, beginning with maybe the most 
relevant for concrete applications, i.e. the multiplication of two delta functions. 

Example nr.l: (S®8)/ a> p) 

First of all we remind that, in this case, the ©, 0^ and Q ex multiplications all coincide. If 
the following limit exists for some a, (3 > 0, we have 



(*©<*)(,*0 (*) = (* = Km 8^(x)8^(x)*(x)dx 



lim n a+l3 / ®(n a x) ^(n 13 x) ^(x)dx, 



for all *$>(x) £ X>(R). It is an easy exercise to check that this limit does not exist, if ^(O) ^ 
and a = (3. Therefore we consider, first of all, the case a > (3. In this case we can write 

(5 5)%»® (*) = n 13 J <&{x) $(ra^ tt ) ^{xn- a ) dx, 

and the existence of its limit can be proved, as in (2J [3], using the LDCT. Choosing §{x) as in 
()3.ip . and defining B m = x m dx, which is surely well defined and strictly positive 

for all fixed even m, it is quite easy to deduce that 

(*0%,»(*)=| o _ a>/»(l + i). <3 ' 12) 

As we can see, this result is quite close to the one in (|2.6p . Analogously to what already stressed 
in Section 2, formula (|3.12p does not imply that (5 8)( a ,/3)(^) does not exist if a < (3 (l + — ) 
because using the LDCT we only find sufficient conditions for (5 <5)( a ,/9)0l0 to exist. However, 
with respect to our results in [2], here we can say a bit more, because we have (6 0i5)( Q m($) = 
((5 <5)(/3,a) 0^)) which was not true in general for the multiplication ®(a,/3) introduced in [2]. 
Therefore we find that 

while nothing can be said in general if a € /3 (l + — ) 1 , /3 (l + 
Example nr. 2: (<5 S')f a m 

In this case the three multiplications 0, 0^ and Q ex do not need to coincide. Indeed we will 
find serious differences between the three, as expected. 
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First of all we concentrate on the computation of (5 Qd8')r a py Because of (j3.9j) this will 
also produce (5 Q ex £')(/3 »)• We have 

(tfOd*')^ (*) =n a+2/3 [ §(n a x)&(nPx)y{x)dx, 

Jm. 

where is the derivative of Again, it is easy to check that the limit of this sequence does 
not exist, if a = (3, for all \&(x) E T>(M) but only for those ^f(x) which go to zero fast enough 
when x — > 0. Let us then consider (5 Qd 8')^'^' (\P) for a > (3. In this case we can write 

(5 Q d 8)<£> P) (tf) = n 2/3 ^ $(z) $'(xn^- a ) $(xn' Q ) dx, 
and by the LDCT we deduce that 

L u ' « ^ p m _ 1 , 

where K m = ^ /_ x x™" 1 $(x) dx. We see that, contrarily to (12. 7D . we can obtain a non trivial 
result with the Qd multiplication. It is clear therefore that also the (80 ex 8') can be non trivial, 
as remarked above . 

The situation is completely different for the (5 5')(/3,a)(^ / )- Indeed, it is not difficult to 
understand that the LDCT cannot be used to prove its existence. The reason is quite general 
and is the following: 

suppose that for two distributions T and S in £>'(R) (T©d5)(a,/3) exists for a and (3 such that 
a > 7/3, where 7 is some constant larger than 1 appearing because of the LDCT. For instance 
here 7 = ^ry, while in Example nr.l we had 7 = 1 + ^- Therefore, since if (S Q ex T)m )0 ,-\ and 
(S Qd T\ a> p) both exist and coincide, using (|3.1ip we have 

(S T) (Q)/3) (*) = \ {(S Q d r) (a ^)(*) + (S Q d T) {M (V)} (3.15) 

Of course, (S TW^^) exists if (S Qd and (5 TWqA^) both exist, which in 

turn implies that a > j(3 and, at the same time, that (3 > 7a. These are clearly incompatible. 
Therefore in order to check whether (S T)( a) m(\I/) exists or not it is impossible to use the 
LDCT which only gives sufficient conditions for the multiplication to be defined: some different 
strategy should be considered. 

Example nr.3: (8' 

As for Example nr.l we remark that here the 0, 0^ and Q ex multiplications all coincide. If 
the following limit exists for some a, f3 > 0, we have 

{8' © 8 f ) {aS (¥) = {8' 8') {M (¥) = J (x) 6™ (x) *(x) dx = 
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lim n 2a+2/3 f &(n a x)&(n /3 x) 9(x)dx, 

rwoo J R 



for all ~$>(x) E D(K). As before, it is quite easy to check that this limit does not exist, if ^(O) ^ 0, 
when a = (3. Therefore we start considering the case a > (3. In this case we can write 



(<*' 5')^ (*) = n a+2/3 J 1 &(x) &(xnP- a ) * 



(m a ) dx, 



and again the existence of its limit can be proved using the LDCT. Choosing <fr(x) as in (|3.ip 
and defining B m = J_j x m ~ l &'(x) dx, which surely exists for all fixed even m, we deduce 
that 



B m 5(V), a = (3 



m+l 



(^OmOH n (3 - 16) 

L U ' ^ m-2' 

However this is not the end of the story, because we still can use the symmetry (S' QS')^ a ^(^f) = 
(5' #'Wa)0^) discussed before. We find 

( 1 0, «</3^T, or a 



while nothing can be said in general if a E 
restrict to the following values of m: m = 4, 6, 8, . . .. 



m- 


fl 


> m 


-2 


m- 


fl 


m4 


1 


m— 


2 



m-2 ' 



. Needless to say, we need here to 



Summarizing we find that results which are very close to those in [2] can be recovered with 
each one of the definitions in (|3,5p . (|3.6p or (|3,7p . The main differences essentially arise from 
the lack of symmetries of these two last definitions compared to definition (|3,5p and the one in 
0. 



4 More spatial dimensions and conclusions 

While the definition of the multiplication given in [2], as we have stressed before, cannot be 
extended easily to M. d , d > 1, definitions (|3.5p . (|3.6p or (|3,7p admit a natural extension to any 
spatial dimensions. We concentrate here only on the symmetric definition, (|3.5p . since it is 
the most relevant one for the application we are interested in here. Of course no particular 
differences appear in the attempt of extending 0^ and Q ex to d > 1. 

The starting point is again a given non negative function E T>(M. d ) with support in 

I\ := [—1,1] x ••• x [—1,1], and such that f, $(x)dx = 1. In this case the delta-sequence 

S v ' 1 

d times 

is dfflix) = re da< 3?(n a x), for any choice of a > 0. The same results listed in the previous 
section again hold in this more general situation. For instance, if T E T>'(M. d ) then {T^ (x) = 
(T * Sn)(x)} is a sequence of C°° functions and it converges to T in D'(M. d ) as n diverges, 
independently of a > 0. 
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Therefore, let us consider two distributions T, S G T>'(R ), and let us consider their convolu- 
tions T^ a \x) = (T*5 { n ] ){x) and si P \x) = (S*8 i J 3) )(x) with 5 { n\x) = n da <£>(n a x), for a, (3 > 0. 
As usual, Tn a \x) and S^\x) are both C°° functions, so that the following integral surely exists: 

(S T)t 0) m = \f \S<?> (x) ltf> (x) + SW (x) rW (£)] *(x) <fc, (4.1) 

V^(x) G As before the two distributions S and T are 0-multiplicable if the following 

limit exists independently of ^(x) G T>(R d ): 

(50T) M) (f)= lim (S0T)W)(f). (4.2) 

We consider in the following the ©-multiplication of two delta functions, considering two 
different choices for the function $(x) both extending the one-dimensional case. 

The starting point of our computation is the usual one: if it exists, (5 6)( a ,f3)( 1 ^) must be 
such that 

(S&S)( a0) (^) = lim n da+d P [ <$>{n a x) $(r/x) f(x) dx. 

n^oo J Rd 

Again, this limit does not exist if a = /3, but for very peculiar functions ^(x). If we consider 
what happens for a > (3 then the limit exists under certain extra conditions. 

For instance, if we take $(x) = nf=i $( x j)i where &{xj) is the one in (13. ip . the computation 
factorizes and the final result, considering also the symmetry of the multiplication, is a simple 
extension of the one in (|3. 13|) : 

[0, a <^(i + i) ora>/3 (i + i.), 

A different choice of &(x), again related to the one in (|3.ip . is the following: 

, v f KTI— eXp{ n .A -i |, ||x|| < 1 , 

$(x) = < ^ ll^ll 1 " (4.4 

1 0, Ibll > 1, 



where N' m is a normalization constant and ||x|| = \l x\ + • • • + x^. With this choice, defining 
C m = -^-^ / R ||x|| <J?(x) and using the symmetry property of 0( aj/ 3), we find 



^ 0, <*<f>{l + £) or Q >/J(l + S. 



Therefore the limit defining the product of two delta can be defined, at least under certain 
conditions, also with this choice of $(x) . The main differences between the above results are 
the values of the constants and the fact that d explicitly appears in the result in (|4.3p , while it 
only appears in the condition relating a and (3 in (|4.5p . 
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The conclusion of this short paper is that the use of sequential completion, properly adapted 
for our interests, is much simpler and convenient. The next step of our analysis will be to use our 
results in applications to three-dimensional engineering structures, trying to extend the results 
in lUE]. 
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